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We study the occurrence of modulational instabilities in lattices with non-local, power-law hop- 
pings and interactions. Choosing as a case study the discrete nonlinear Schrodinger equation, we 
consider one-dimensional chains with power-law decaying interactions (with exponent a) and hop- 
pings (with exponent /3): an extensive energy is obtained for a,f] > 1. We show that the effect of 
power-law interactions is that of shifting the onset of the modulational instabilities region for a > 1: 
at a critical value of the interaction strength, the modulational stable region shrinks to zero. Sim- 
ilar results are found for effectively short-range non-local hoppings (/3 > 2): at variance, for longer 
ranged hoppings (1 < /3 < 2) there is no longer any modulational stability. We also discuss the 
stability regions in presence of the interplay between competing interactions - e.g., attractive local 
and repulsive non-local interactions. We find that non-competing (competing) non-local interactions 
give rise to a modulational instability emerging for a perturbing wavevector q — tt (0 < q < n). 
The hopping instability instead arises for g = perturbations, thus the system is most sensitive to 
perturbations of the order of the system's size. Since for a > 1 and /3 > 2 these effects are similar to 
the effect produced on the stability phase diagram by a finite-range interactions and/or hoppings, 
we conclude that the modulational instability is "genuinely" long-ranged for f < /3 < 2 non-local 
hoppings. 

I. INTRODUCTION 

The investigation of the effects of the interplay between discreteness and nonlinearity is a long-standing 
argument of research in the study of the dynamical properties of nonlinear lattice models [IHTj . A typical 
feature exhibited by nonlinear classical Hamiltonian lattices is the existence of discrete breathers, i.e. 
time-periodic and localized in space solutions of the equation of motions. The study of their dynamical 
stability, as well as their robustness in long transient processes and thermal equilibrium, has been the 
subject of an intense experimental and theoretical work [7]. 

The interplay between discreteness and nonlinearity is also crucial for the occurrence of modulational 
instabilities (MI), well known in the theory of nonlinear media [1] [2]. MI are dynamical instabilities 
characterized by an exponential growth of arbitrarily small fluctuations resulting from the combined effect 
of dispersion and nonlinearity. The occurrence of modulational instabilities has been studied in a number 
of physical systems, ranging from fluid dynamics [S] to nonlinear optics [9j. The role and the consequences 
of the MI in the dynamics of discrete systems have been extensively studied: the MI was discussed in 
the context of the discrete nonlinear Schrodinger equation (DNLSE) [10 , which is a paradigmatic lattice 
model used to study nonlinear discrete dynamics [31 S] . The DNLSE is commonly used to describe 
the effective dynamics in different physical systems of interest, including the dynamics of ultracold atoms 
in optical lattices [12] and optical waveguide arrays [T^ . For ultracold bosons in optical lattices the onset 
of MI was analytically predicted [14) and experimentally observed (15j , and in nonlinear waveguide arrays 
the experimental observation of the MI was also reported [T5] . 

In this paper we study the occurrence of MI in the DNLSE in the presence of non-local long-range 
hoppings and interactions. A motivation for such a study comes from experiments with ultracold dipolar 
bosonic gases [T71 [TS] which have been Bose condensed in the last years by several groups [THHH] , from 
the attainment of quantum degeneration for ensembles of polar molecules [22] . and from the recent 
experimental investigations of strongly interacting Rydberg gases [23H25[ . Since the interaction potential 
between (di)polar atoms (molecules) decay as a power law 1/r^ (for Rydberg gases interacting through 
Van der Waals interactions as 1/r^), recent experiments with dipolar gases in optical lattices |26j, as well 
as the realization of long-lived dipolar molecules in a 3D periodic potential [37] and in perspective also 
the dynamics of Rydberg atoms in optical lattices [28l [29] , open the possibility to study DNLSE with 
non-local interactions. 

Our other motivation is related to the wide interest for systems with long-range interactions |31| : 
in these systems the range of interaction of the constitutive units is not bounded. A typical form of 
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interactions, relevant for a number systems ranging from gravitational ones to dipolar magnets and 
gases, is provided by the power-law decay Xjr'^ (e.g. for gravitational systems 7 = 1) where r is the 
distance among the constituents. For statistical mechanics models, like the Ising or more generally the 
0{n) models, the possibility to have power-law couplings make possible the appearance of a rich phase 
diagram [3ll4]. 

A first criterion to determine the long-rangedness of a system with power-law decaying interaction 
is the comparison with the dimension d of the space; as 7 is smaller or equal than d, if the system is 
homogeneous and the interaction favours homogeneity we obtain a diverging energy density, thus in order 
to obtain a well defined thermodynamic limit (if relevant) a rescaling of the energy is in order (the so 
called Kac rescaling) [35 . In the following we will refer to this region as the non-extensive long-range 
region. If 7 is larger than d the energy of the system is extensive and it is usual to individuate a value 
of 7, which we denote by 7*, such that for 7 > 7* the system behaves as a short-range system. Since it 
is 7* > d, there is a region of values of 7, given by d < 7 < 7*, in which the behaviour of the system 
significantly differs from the properties of the same system with short-range interactions, although the 
energy is extensive. Such region is the extensive long-range region, also referred to as the weak-long-range 
region [35], where "weak" refers to the extensivity of the energy. The actual value of 7* depends on the 
specific model and the dimension: for 0(n) models in d = 1 it is 7* = 2 [551 f57| . 

Both the thermodynamics and the dynamics of long-range interacting systems are extremely interesting 
[311135]: in particular, in the long-range region the dynamical evolution evidences that the system may stay 
in a quasi-stationary metastable state (different from the thermal equilibrium one) for a time exponentially 
growing with the size of the system. Such metastable state is reached after a short-time dynamics, referred 
as violent relaxation [35] . 

While the main body of nonlinear wave systems studies have dealt with short-ranged systems, some 
extensions of these results for long-ranged systems have appeared in the literature. Among theoretical 
studies addressing the properties of discrete systems with different kinds of nonlocal dispersion or nonlocal 
nonlinear interaction we mention |38H42) which focus on existence and stability of localized excitations. 
For a recent review on these and other problems emerging in nonlinear lattice systems theory see |30j . 

The purpose of the present paper is to clarify how modulational instabilities emerge in nonlinear lat- 
tices with non-local interactions and hoppings, aiming both at unveiling if (and in what conditions) short 
time dynamical instabilities occur in nonlinear lattices and at clarifying the nature of the emerging mod- 
ulational instability. We choose the DNLSE as a case study not only due to its paradigmatic usefulness, 
but also due to its relation with XY [i.e., 0(2)] models: when number fluctuations are frozen, the kinetic 
term in the DNLSE energy is basically the XY model (see the discussion in Section |ll]). This is the reason 
why we choose to consider not only power-law interactions (as it is relevant for experiments with ultra- 
cold dipolar bosons in optical lattices), but also power-law hoppings (which corresponds to power-law 
couplings in 0(2) models). Using the DNLSE we study the modulationally stable and unstable regions 
in presence of power-law non-local interactions and hoppings, discussing also the interplay between local 
and non-local interactions, e.g. local attraction and non-local repulsion. 

The plan of the paper is the following: in Section |ll| we introduce the DNLSE with long-range hoppings 
and interactions and discuss its relation with other statistical mechanics models. In Section [IIII we derive 
the Bogoliubov spectrum of elementary excitations, present the general framework for the determina- 
tion of the stability regions in presence of power-law interactions and hoppings and specialize it to the 
analysis to the short-range limit reminding the known results of the MI analysis [lOl E] . Our findings 
for power-law interactions are presented in Section [IV| where we also consider the case of attractive and 
competing local and non-local interactions. Sectionfv) deals with a system with non-local hopping and 
local interaction which presents some peculiar features with respect to the long-range interaction which 
are further investigated in |VI[ Our conclusions are in Section [Vll[ 



II. THE DNLSE WITH LONG-RANGE INTERACTIONS AND HOPPINGS 



The DNLSE with non-local interactions and hoppings reads: 

m m 

In Eq. ([T]) T is the time and the indices j, m denote the sites of a lattice. For simplicity we assume that the 
lattice is one-dimensional, but the subsequent analysis can be extended to higher dimensional lattices. 
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The indices j, m then assumes the values j,m = 0, ■ ■ ■ ,L — 1 (L is the number of the sites, taken to be 
even) : periodic boundary conditions wih be also assumed, so that the wavefunction satisfies the condition 
= ijj+L- The Hamiltonian corresponding to ([T]) reads 

^DNLSE = - X! ^pj,mi^in +\^ \ V'j ? ^],m I V'm P ■ (2) 

We denote the diagonal interaction by 
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j.m j,m 



and the next-neighbour interaction as 



Vj,i = U, (3) 



V,.,+i = V. (4) 



The interaction coefficients V^,m iu ([T]) are assumed to be power-law decaying with exponent a, i.e. 

1/ I m ^ i I". Since Kjij = ^,m: implementing the periodic boundary conditions amounts to require 
that Vj^jn = V(j+n) mod L,(m+n) mod L, whcrc j < TU and n = 1, • • • , L - 1. Wc have therefore 

(U m = 0, 

Vo,nr^\^ m=l,...,f (5) 
I (L^ m=f + 1,...,L-1 

The non-local hopping rates tj_rn will be also assumed to be power-law decaying with exponent /3: we 
consider vanishing diagonal hopping {tj,j = 0) and a nearest-neighbour hopping 

= (6) 
Therefore, with periodic boundary conditions we have 

m = 0, 

to.rn = { - ^1 (7) 



m==4f + l,-- - ,L-1 



(L-m)f ~ 2 

Since we want a finite expression for the energy per particle, we will consider 

a>l;/3>l. (8) 

To treat the cases a < 1 or /3 < 1 one should do a Kac rescaling, as it is usually done in statistical 
mechanics models with non-extensive long-range interactions [35] : e.g., for /? < 1 one has to perform the 
substitution 

tj,m ~> ]— • (9) 

A non-extensive ground-state energy is found in our case by a — ^ 1 and/or /3 — >■ 1: the region 1 < /3 < 2 
is the weak-long-range hopping region. We remind that such region is of high interest in the study of 
statistical mechanics models with long-range couplings. Let us consider an Ising model of the form 

-falsing = - ^ JijSiSj {Si = ±1) (10) 
id 

with power-law couplings Jij oc 1/ \ i ~ j p . As usual, i,j denote the sites of a lattice with dimension d. 
It is well known that for nearest-neighbour couplings (formally corresponding to 7 = 00), then there is 
order at finite temperature only if d > 2 However, if the interactions are sufficiently long-ranged it 
is possible to have a phase transition at finite temperature between a ferromagnetic and a paramagnetic 
phase [44| , l45 ] even for d = 1: if 7 > 2, then the critical temperature Tc is vanishing (as in any short-range 
Ising chain [13]), while for 7 < 1 the energy is non-extensive. After the Kac rescaling one easily sees that 
for 7 < 1 there is a phase transition having the critical exponents of the mean-field universality class (see 
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e.g. [46]). For 7 between 1 and 2 (weak-long-range region) there is a phase transition at a finite critical 
temperature. For 7 = 2 a Kosterlitz-Thouless transition takes place [45^ 47^ 48, 

Before moving on to the derivation of the Bogoliubov spectrum of elementary excitations, we pause here 
to discuss the relation between the DNLSE ([T]) and a model widely used in the treatment of long-range 
systems, i.e. the Hamiltonian Mean Field (HMF) model We observe that the DNLSE Hamiltonian 
Q can be written as the sum of a kinetic term -fffci„ and an interaction term Hint- The kinetic part reads 
Hkin — — m '4'j'tj,m'4'm'- by Writing ipj in terms of the local density rij and phase 9j, i.e. tpj — y/nje'^^^ 
the kinetic part reads then H^in = — J2j m ^j,m\/'''''j'^m cos {6j — 0„i). One sees than when the number 
fluctuations are frozen, i.e. Uj ~ n (where n is the average number), the DNLSE Hamiltonian reduces to 
the potential energy (and it has the same equilibrium properties) of the HMF model |35j 

-^HMF = — ^ Jj,7n COS {9j — 9m), 

(where Jj^m = tj^mn), which is nothing but a long-range XY model. This result is of course expected in 
the sense that interacting bosons on a lattice are in the XY universality class, and interacting bosons 
with long-range hoppings have to be in the universality class of the long-range XY model. 



III. MODULATIONAL STABILITY ANALYSIS IN PRESENCE OF LONG-RANGE 

INTERACTIONS AND HOPPINGS 



We study in this Section the Bogoliubov spectrum of elementary excitations, describing the energy of 
small perturbations with (quasi)momentum q on top of a planewave state with (quasi) momentum k |10j . 
The final part of the Section is devoted to briefly remind the results of the short-range limit {V — and 
j3 — >■ 00, i.e. only nearest-neighbour hopping). 

The stationary solutions of Eq. ([T]) are planewaves 

V'j(t) = V'o exp [i{kj - i^t)] : 
V is the chemical potential given (for L oo) by 

hv = -2Up{k) + p[U + 2VC{a)]. (11) 



In Eq. (Ill p is the planewave density {p = \ tpo I ); furthermore (a) is the Riemann zeta function 



oo 



(12) 



and we introduced the function 



e^ik) = J2 



cos (mfc) 



(13) 



Useful properties of the function (13) are recalled and discussed in Appendix [A| 

The stability analysis of planewaves stationary solutions can be carried out by perturbing the carrier 
waves as 

^j{t) = [V^o + u(T)e'« + w*(r)e-*«] e^C^J-"^*). 
Retaining only terms proportional to u/ipo and w/V'o, one gets 



ih 



dr \ V 



A C 

-C* -B 



(14) 



The quantities A, B and C in Eq. ( 14 ) are defined by 



A^2t [Ip {k)~£p {k + q)]+pVg, 



(15) 



5 



B^2t [if, (k) -£p{k^ q)] + pVg (16) 

and 

C = ^o'n- (17) 
In the previous expressions, Vq is the Fourier transform of the interaction ([s]): for finite L it is 



m=0 

For L oo one obtains 



Vq = U + 2V£^ {q) . (18) 



From Eq. (14 1 it follows that the excitation spectrum (i.e., the Bogoliubov dispersion relation) for the 



DNLSE with power-law hoppings and interactions is 

huj± = ^ ^ ± 2tVi: (19) 

where u!± = i^±{k; q) is a function of k and q (respectively momentum of the perturbed and perturbing 
planewaves) . Furthermore 

T = i {2£p{k) -i^ik + q)- £p{k - q)f + 
+ ^ {2£p{k) - £p{k + q)- £p{k - q)} = 

^T{k-q)I^F{k-q) + P^'^ (20) 

where we introduced the function J-{k\ q) = {2i?/3(fc) — ^^(A: + q) — ^/3(fc — q)} /2. In the following we will 
use the convenient dimensionless parameters 

U^^;V^^. (21) 
t t ^ ' 



In terms of the parameters J7, V, the quantity pVq/t entering (201 reads 

PVq 



^ U + 2V£a.{q). (22) 



The carrier wave becomes modulationally unstable when the eigenfrequencies uj± in Eq. ( 19 1 acquire 
a finite imaginary: the condition for stability is then I > 0. A momentum k is then modulationally stable 
if for each q the eigenfrequencies uj± are real, otherwise if it exists a q such that I < then k will be 
modulationally unstable. Since the eigenvalues are unaffected by substituting k with —k and q with —q, 
we will consider k and q both belonging to the interval [0, tt]. Notice that for 5 = it is I = uj± — 0. 

When a momentum k is modulationally unstable, there will be some momentum q for which w+(fc; q) 
and id-{k; q) have an imaginary part '^uj±{k; q). For those values of k and q we write 

r(fc;q)=|3c^±(fc;g) I (23) 

to quantify how much unstable is the perturbed planewave (notice that the imaginary parts of u)j^{k]q) 
and uj^(k]q) are by definition opposite in sign and equal in modulus). For k unstable we will use the 
notation Tj^axi^) — maxg r(fc; q), where the max is taken on all the q such that uP'{k; q) < 0. The value of 
q for which the maximum value of r(A:; q) is obtained will be denoted by qmax (i-e., Tmaxik) = T{k; qmax))- 
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FIG. 1: Stable (white) and unstable (in color) regions for the short-range DNLSE {V — and only nearest- 
neighbour hoppings) for U = 2.5 and U — 1.5 (panel a and b respectively) in the k (carrier wave momentum) 
q (perturbing wave momentum) plane. The absolute value of the imaginary part of the frequencies r(fc; q) is 
depicted (from light to dark red as T{k;q) increases). The continuous lines are equispaced isolines whose spacing 
is set to 0.2. The thick red line indicates the position of qmax- 



A. The short-range limit 

In this Section we review the results and the region of stability for the short-range limit, having only 
local interaction {V ~ 0) and nearest-neighbour hopping: <jj±i = t and tj^„i = for m =/: j ± 1 [formally 
this is the hmit ^ — > oo in Eq. Q]. 

It has been shown in [TU] that the onset of MI occurs at kcr — tt/2: i.e., the momenta k < tt/2 are 



modulationally stable, while for k > 7r/2 are unstable. The quantity I, defined in (20) and giving the 
stability regions reads 

T = 4cos2fcsin''--|-2t7cosfcsin2- : (24) 

2 2 ^ ' 

one readily sees that the momenta k < tt/2 are stable. For k > 7r/2 all the momenta k are unstable, 
irrespectively from U. However, one sees that for k > tt/2 and U > 2 each q is unstable, while for [/ < 2 
there are stability regions: these stable regions are found to be bounded by the regions in which q is 
between qcr and tt and k is between kcr and tt, where 

I. ^ 
krr = TT — arccos — 

2 

and 

qcr = 2 arcsin 

The resulting plots of stable and unstable regions for U > 2 and U < 2 are in Fig. [l] where we plot as 
well as contour plot the value of T{k;q) in the unstable regions (the larger F, the more unstable is the 
dynamics). In Fig. [2] we also plot as a solid line the values qmax{k) for which the maximum value of F, 
at the fixed k, is reached. 



IV. POWER-LAW INTERACTIONS 



In this Section we consider the case of power-law interactions in presence of (local) nearest-neighbour 
hoppings {tj.j±i — t and = for m 7^ j ± 1). One has then 

1 = 4 cos^ k sin'* | -I- 2 cos fc sin^ | (t^ + "iVlo, (<?)) ■ (25) 

We consider only the cases where V > and the local interaction U can be either or positive or negative. 
The cases where V is negative can be easily derived from the cases where V > Ohy noticing that under 
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FIG. 2: Stability in the k-q plane (see caption of Fig. [T| for power-law interactions in presence of local nearest- 
neighbour hoppings with a = 1.5 and U = 2.5. We plot T{k;q) for the four values of = 0.2, 0.5, 2 and 4 in 
panels a, b, c and d respectively. 



the following transformation: 



cos(fc) 
U 
V 



cos(fc) 
-V 



(26) 



one obtains the same dependence on q and a of the stability conditions that one has when V is positive. 
E.g., if for a fixed value V = Vq > and U = Ua < one has that the momentum k is stable, then 
momentum tt — k will also be stable for V — —Vq < and U = —Ua > 0. 



A. Repulsive interactions: U > 0, V > 

We consider here U and V to be positive. Since d£a/dq < for a > 1 and q G [0, tt], one can shows 
that: 

• for 

tJ > 2VC{a){l - 2^-") 

the momenta k < tt/2 are stable; 

• for 

U < 2VCia){l - 2^-") - 2 

the momenta k < tt/2 are unstable. 
It follows that the critical value k^r as a function of V is given by 

U 



cosfccr = ^C(a)(l - 2'"") - 



(27) 
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FIG. 3: kcr vs V from ( |27[ ) (in the figure U = 1 and a = 3). The dots represent values obtained by direct simulation 
of the DNLSE, see the text for details (errors are smaller than the symbols). 



i.e., for 2t^C(a)(l — 2^^") — 2 > C7 one has kcr — 0. For k > tt/2 it is easy to see that all momenta are 
unstable against perturbations at q — 0. 

The instability regions are depicted in Fig. [2] as V is increased - for C7 > 2, one can identify four regions: 
for 2V I £a{n) | +2 < C7 with £„(7r) = C(a)(l - 2^-"), then the moment,a k larger than tt/2 (smaller 
than 7r/2) are unstable for all q. Increasing V, one has that for U — 2 < 2V \ laiT^) \< U & stable region 
forms around k = n,q = n, while for U < 2V | £Q(7r) |< [7 + 2 an unstable region appears close g = tt for 
momenta k between kcr given by Eq. (27) and 7r/2. When U + 2 > 2V \ iaiT^) I) then all the k becomes 
unstable and the instability starts from q = tt. Therefore, the planewave kcr is rendered unstable by the 
wavevector q = tt, thus the system is the most sensitive to short wavelength perturbations. 

The behaviour of fccr is plotted in Fig. |3] the analytical prediction (27 1 (valid for a > 1) is compared 
against numerical findings obtained numerically solving the DNLSE. The numerical solution has been 
obtained on a finite size system (we choose L = 512) whose coherence has been monitored by inspecting 
the absolute value of the following order parameter [TJ] 



(28) 



where tpk{T) — '4'in{T)e is the Fourier transform of the wavefunction. The initial wavefunc- 

tion ipjir = 0) is chosen as a planewave with wavevector k perturbed by the highest frequency wavevector: 
'(/'^(t = 0) = e^^^ -\- ee"^^ , with q = tt (notice that qmax = for non-competing interactions, as one can see 
from Fig. [2]) . The ratio between the amplitudes of perturbing and perturbed wavefunction is set to be 
e = 10"'^ with p = 1. As we can see in the left panel of Fig. Wl if prepared in a modulationally unstable 
initial state, the system loses coherence after a time which diverges as we approach the momentum kcr'- 
denoting with t/ the time after which the instability is observed using the order parameter ( 28 1 , and 
noticing that the quantity (20) is vanishing (at q — tt) as ^ {k — kcr), one can estimate kcr from the 
numerical data using the dependence t/ oc {k — kcr) (see Fig. |4j right panel). This divergence has 
been used to extract the numerical values of fccr shown in Fig. [3] 
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FIG. 4: Time evolution (on the left) of the modulus of the order parameter </?, defined in equation (28 1, for five 
different values (from left to right k = 77, 76, 75, 74, 73, 72 in units of 27r/L with L = 512) of carrier planewave 
wavevector for the DNLSE with parameters U — 1, V = 1.25 and a — 3. The lifetimes of the initial states are 
depicted in the right part of the figure. These times have been fitted with a function r = const. * (k — kc 
(dotted line) to obtain the a numerical estimates of k^r (vertical line). 
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B. Competing interactions: U < Q, V > Q 

For fc < ^ a necessary condition for stability is given by |C/| < 2V(^{a) that is obtained by analyzing 
the stability at g = 0. The critical momentum is given by 



mm 

ge[0,7r] 



arccos ^7— — 

V 2sin2(|) 



(29) 



all momenta k < k^r are stable. 

An important feature of the case with competing interactions is that the most unstable perturbations 
can arise for a value of q* different from and tt (i.e., < g* < tt). In the following we determine for 
what conditions at the critical value kcr there is an instability at a qmax = 9* 7^ 0, tt. We will refer to 
these values of q* as to finite values for the occurrence of modulational instability: this is because the 
instability develops on a length scale ~ l/?*. It is intended that if the instability arises at g = tt, this 
develops on a length scale of the order the lattice unit, while an instability at g = involves length of the 
order of the lattice size: the latter will be the case for non-local hoppings with weak-long-range exponents 
1< ^ < 2. 

An example of a finite q* is shown in Fig. [5] To see how this can arise we analyze for simplicity the 
situation at fc = 0, which can be easily generalized to a finite value in the interval < fc < 7r/2. In order 
to observe an instability region like the one in Fig. [5] we have to impose that the equation 

h{q) = \U\ (30) 

with 

%) = 2sin(|)'+21/£„(g), (31) 



obtained from Eq. (25) by the substitution k = Q, has exactly one finite solution for fixed values of the 
parameters (a, ?7, yjand no solutions for smaller values of tj . After some algebra one can derive a set of 
conditions on the coefficients (a, V) such that the curve h[q) possesses one minimum. It turns out that, 
if J7 = minqg(o,7r) h{q)^ then there is a S (0, tt) such that the instability region is tangent to the fc = 
axis. For smaller values of U the previous equation h{q) = \U\ does not have any solution and one obtains 



a finite value of kcr given by Eq. ( 29 1 



We can thus assert that the presence of competing interactions may originate the wavelength q* smaller 
than TT (and larger than tt), unlike the case of non-competing interaction examined in the previous 
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FIG. 5: Stablity regions (see caption of Fig. [T]) for the case of competing interactions with a = 2, U — —0.7 and 
V = 0.5. The value of the most unstable q which defines kcr is denoted by q* . 

subsection for which the system is most sensible to perturbations at g = tt. This is a general feature of 
systems with competing interaction acting on different scales which, if properly tuned, give rise to the birth 
of a new intermediate lengthscale (in this case of the order of l/q*)- For similar phenomena ultimately 
leading to stripe formation and more generally spatially modulated patterns in different contexts see 
e.g. Hi [SO]. 

For completeness we list the conditions on the parameters (a, V) such that the function h{q) has a 
minimum for a finite value q* of the perturbing wavevector: 



for 1 < a < 



for a* < a < 3 



for a > 3 



^ < (2-2i--)C(.) 



2(l-23-")C(a-2) ^^^^ 



^ < < W7 ZT. (34) 



2|C(a-2)| 2(l-23-")C(a-2)- 



(similar results are found for < fc < 7r/2). The value a* in (33) is given by « 1.513, the unique solution 
of the equation: 

2(1 - 2'-'^)aa - 2) - (1 - 2i-")C(a). (35) 
The condition ensuring the stability for every k £ [§, tt] is given by: 

\U\ > 2VC{a). (36) 
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FIG. 6: Stability regions (see caption of Fig. [TJ for the case of long-range hopping and local interaction U = 1 
and = 0. The panels a, b, c and d refer to the values of /? = 4, 3, 2, 1.5. 



V. POWER-LAW HOPPINGS 

In this Section we consider the case of power-law nearest-neighbour hoppings, with exponent j3 > 1: 
one has 

I=\ {2lp{k) - ep{k + q)- £p{k - q)f + 
+ I {2£p{k) - ip{k + q)- ip{k - q)} = 

^F{k-q){F{k-q) + U) (37) 



While in the cases considered in Section IV the instabilities arise in the higher frequency range oiq[q = tt] 
for non-competing interaction or at a finite value oi q[q = q* G (0, 7r/2)] for competing interactions, with 
non-local hoppings even the long-wavelength perturbations can affect significantly the stability properties 
of the system. This can be verified by inspection of the behaviour of I for small q: one finds for g — ^ 

Uq^d^ip 
4 9fc2 ■ 

The investigation of the behaviour of the second derivative of tp{k) reveals that d'^lp/dk'^ is positive for 
1 < f3 < 2 for each k. It follows that 

fccr = for 1 < /3 < 2, (38) 

i.e. the modulational stability regions shrinks to zero in the weak-long-range regime, irrespectively from 
U and V (we assume for simplicity in this Section U > and V > 0). The same analysis shows that for 
y = the critical value kcr does not depend on the specific value oi U ior (3 > 2. The above scenario is 
confirmed by Figure |6] where we plot r{k;q) as /3 is increased. The values of kcr as a function of /3 are 
plotted in Fig. [t] Notice that kcr tends to 7r/2 when the hopping exponent approaches the short range 
limit /3 ^ oo. As done in Section [TV] we compare the analytical results with numerical simulations of the 
DNLSE obtaining a very good agreement. 
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FIG. 7: Critical momentum fee 
numerical data. 



vs the power-law exponent /3 of power-law hoppings with = 0. Crosses are 



It should be stressed that the hfetimes of the modulationally unstable states in this case is typically 
longer than the ones encountered in Section |IV[ This will be supported in the next Section where we 
directly compare instabilities arising from the long-range interaction and hopping respectively. 

When we introduce the long-range interaction V ^ the instability due to long-range hopping remains 
unchanged since it is due to the vanishing of q) (see Equation ([37])), while the instabilities already 
discussed in Section [TV] are possibly generated. In this setting the contribution of both instabilities has to 
be taken into account, with g = instabilities generated by long-range hopping and q = n (or q = q* < tt 
if the competition is sufhciently strong) instabilities due to interaction. We comment on these scenarios 
in the next Section. 



VI. COMPARISON OF INSTABILITIES ARISING WITH NON-LOCAL INTERACTIONS 

AND HOPPINGS 



The hopping instability described in Section |V] exhibits an important difference with the one arising 
from interaction described in Section IV since kcr becomes unstable for q — perturbations, i.e. with a 
size of the order of the system's size. 

In order to compare the timescales on which these two kind of instabilities act we have considered 
a case where we switch on and off alternatively the long-range interaction and hopping (see Figure [s]). 
We prepare the two systems with a planewave k slightly inside the instability region perturbed by a 
q equaling tt and 27r/L (the widest available perturbation with a nonvanishing imaginary Bogoliubov 
frequency) . By inspecting the contour plots (panel a and b) we foresee a longer lifetime for the instability 
induced by the non-local hoppings: this result is confirmed by numerical simulations, shown in panel c 
of Figure |8] This is a general feature of the long-range hopping instability: we observe that at the the 
critical value k^r, by definition, T{kcr',qmax) = both for non-local interaction and hoppings. However, 
entering the unstable region gives a vanishing value of T{k; q) for g = and k > kcr and r(/e; q) oc q for q 
small for long-range hoppings. At variance for long-range interactions as k is slightly larger than kcr then 
r(fc;'7mst) generally acquires a finite value at the wavevector qinst at which the instability arises (with 
Qinst = for the non-competing case and qinst = 9* G (0j''''/2) for the competing one). 

Let us examine the peculiarities of the instabilities due to interaction and hopping and examine whether 
these may be reproduced with finite-range couplings or hoppings. As for long-range interaction is con- 



cerned, in the non competing case (examined in Subsection IV A I , we observe that the situation is not 
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FIG. 8; Panel a and 6 are contour plots of r(/c; q) (see caption of Fig. [T| for U = V — 1 and a — 2.5, jS — oo 
(long-range interaction) and a = oo, /3 = 2.5 (long-range hopping). In panel c we depict the time evolution of 
the order parameter defined in Eq. (281 for the initial conditions indicated in the contour plots with points. As 



we can see the long-range hopping instability (green line) takes a much longer time to set in than the long-range 
interaction instability (green line). 



a) n 




b) 




1 

0.8 
0.6 
0.4 
0.2 




FIG. 9: Contour plots of r(fc; q) (see caption of Fig.[T]) for models including only on-site interactions ([/) and next- 
and next-to-nearest-neighbour interaction {V and Vjj+2 = V2) and hoppings {t = 1 and = ^2) respectively. 

Panel a refers to values U = —0.8, V — 0.4, V2 = 0.05 and t2 = while panel 6 refers to values U = 1, V — V2 ~ 
and t2 = 0.2. 



very different from the finite-range one with nearest- neighbour interaction (a — >■ cxi): we can indeed find 



an effective nearest-neighbour V interaction giving rise to the same kcr (see Equation (27)) 



V ^VC{a){l-2^'°'). (39) 



Similarly in the case where competition in present (described in Subsection IV B ) the case with a finite 
can be seen to be similar to the one obtained with a 00: in fact one can generate a new wavevector 
q* even when a — >■ 00, and thus the long-rangedness of the interaction is not actually playing a major 
role. As an example in Figure [9] (panel a) we have considered a case where a nearest-neighbour and 
next-to-nearest-neighbour interaction generates a stability diagram highly resembling the one shown in 
Figure [5j 

When we move to the hopping instabilities (examined in Section |v| we observe that while a finite-range 
hopping can generate q = instabilities for /c < tt, the value of k^r is effectively limited by the range of 
the hopping and it is not 0. Indeed one can convince himself that for a hopping of range R (i.e. tij = 
for \ i — j \> R) then the critical momentum kcr can become small as tt/ {2R) due to the q — instability. 
Thus the case 1 < /? < 2 where kcr = is singled out as a case where the instability is genuinely due 
to long-range nature of the hoppings. In panel b of Figure [9] we depict the stability diagram of a case 
with nearest- and next-to-nearest-neighbour hopping mimicking the effect of long-range hopping shown 
in panels a, 6 of Figure [6j with a finite value of kcr- 
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VII. CONCLUSIONS 



In this paper we studied the occurrence of modulational instabihties in nonlinear lattices with long- 
range hoppings and interactions. We were motivated by experiments of (di)polar gases in optical lattices 
and by the interest in the study of dynamical regimes in systems with long-range interactions. Using the 
discrete nonlinear Schrodinger equation in one dimension, we considered power-law decaying interactions 
(with exponent a) and hoppings (with exponent /3). We showed that the effect of long-range interactions 
is that of shifting the onset of the modulational instability region for a > 1 (corresponding to an extensive 
energy): at a critical value of the interaction strength, the modulational stable region shrinks to zero. 
Similar results are found for short-range non-local hoppings (/3 > 2): at variance, for longer-ranged 
hoppings (1 < /? < 2) there is no longer any modulational stability. 

Instabilities due to interaction generally differ from the ones due to hopping since the first ones, are 
sensitive to finite wavelength perturbations while the latter are most sensitive to perturbation of the 
order of the system size. Such hopping generated instability turns out to have generally longer timescales 
than the interaction generated ones. If we allow interactions acting on different scales to compete we may 
generate instabilities with longer, but finite, wavelengths, in analogy with what is met in other systems 
with competing interactions [3T] . 

The instabilities met in the long-range interacting and long-range hopping for /3 > 2 are not specific to 
the long-ranged nature of the interaction or hopping; in fact their effects are in principle not different from 
finite range cases. As for very long-ranged (1 < /? < 2) hopping is concerned, we found that it gives rise 
to genuinely long-range instabilities, since they cannot be reproduced with suitably chosen finite-range 
hoppings. 

Acknowledgements We thank M. lazzi, A. Smerzi, G. De Ninno and F. Staniscia for very useful discus- 
sions. 



Appendix A: Useful properties of la{k) 



The analysis of stability regions presented in the main text is based on the study of the quantity I 
defined in (20), which in turn contains the function ^a{<l) defined in Eq.(13|: 



COS {mq) 



(Al) 



(with a > 1). We are interested in the domain q € [0, vr]. From the definition it follows that £a{0) = C{ct)', 
it is also 



C(«). 



(A2) 



The plot of £aiO) and ^a(7r) as a function of a is drawn in Fig. [loj notice that the behaviour of £a{0) for 
a — >■ 1 and a — > oo is given respectively by limQ_i.i £^(0) — oo and limQ,_j.oo ^a(O) = 1. For iaij^) one has 
limQ,_>.i ^c((7r) = ln2 < and Imia^ao (■a{'^) = ^1- 

The derivative of has a different behaviour forl<a<2,a = 2 and a > 2. It is 



(^) = 



for a > 1 and 



9^ 
dq 



(0) = 



a > 2 

-f « = 2 
-oo 1 < a < 2 



(A3) 



The second derivative of i(a) can be computed explicitly and it gives: 



dq^ 



{q) = -lc.-2{q). 



(A4) 
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FIG. 10: Plot ol £a{0) (solid line) and lain) (dashed line) as a function of a. 



For 1 < a < 2 then (q) is a positive function and we have (q) = ^, constant over q e [0, tt]. The 
second derivative takes the following values at the extrema of the interval [0, tt] for a > 1: 

{oo a <2 

1/2 « = 2 

-oo 2 < a < 3 ^^^> 
-({a -2) a > 3. 

Finally at g = tt we have (tt) = (1 — 2^^°')({a) for every a > 1. 
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